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Abstract 

We construct the covariant effective field theory of gravity as an expansion in inverse 
powers of the Planck mass, identifying the leading and next-to-leading quantum corrections. 
We determine the form of the effective action for the cases of pure gravity with cosmological 
constant as well as gravity coupled to matter. By means of heat kernel methods we renormalize 
and compute the leading quantum corrections to quadratic order in a curvature expansion. 
The final effective action in our covariant formalism is generally non-local and can be readily 
used to understand the phenomenology on different spacetimes. In particular, we point out 
that on curved backgrounds the observable leading quantum gravitational effects are less 
suppressed than on Minkowski spacetime. 


1 Introduction 

In recent years the fact that quantum gravity can be treated as an effective field theory (EFT) 
has become an established fact [H El El n [5]. The modern perspective is that even the stan¬ 
dard model of particle physics, originally constructed under the guiding principle of perturbative 
renormalizability, must be seen as an EFT [HlElEli. The same applies to QED, theory that 
claims the best comparison between theory and experiment, together with many other success¬ 
ful quantum field theories. Thus, in complete opposition to what has been thought for a long 
time, gravity can be treated as a perturbative (effective) quantum field theory alongside with the 
other fundamental forces. It is now also recognized that there is no fundamental inconsistency 
between general relativity and quantum mechanics, at least at low energies: the EFT of gravity 
is a perfectly well defined and predictive quantum theory of gravity. While the standard model 
breaks down at a given UV scale, the EFT of gravity will break down at the Planck scale. Any 
UV completion of gravity will have to reproduce the EFT predictions at sufficiently low ener¬ 
gies. This is one reason why we should work out the EFT of gravity. On the other side, as the 
example of QCD has shown clearly, even if we know the UV completion of a theory, so that all 
physical questions are mathematically well posed, it can still be extremely difficult to work out 
the low energy spectrum and the effective action. Thus even if we knew today the fundamental 
description of gravity at the Planck scale, we will probably still need to resort to EFT techniques 
to work out the accessible phenomenology. 

To date, the EFT of gravity has been developed on a Minkowski background, using Feynman 
diagram techniques laiiniiiiiiis], or their modern versions [lailllllsllle], and power counting 
arguments HZj. The paramount result is the calculation of the leading quantum gravity cor¬ 
rections to the Newtonian potential [1]. This is probably the most reliable result we have in 
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quantum gravity so far, but these corrections are very small and a possible direct observation of 
their effects is for the moment quite remote. This fact in turn explains why Einstein’s theory 
of general relativity (GR) is so successful, and predicts that it will still be for a broad range of 
energy scales. The smallness of the quantum effects is the result of the wide separation between 
the scales where observations are made and the Planck scale. This has lead to state, somehow 
paradoxically, that the EFT of gravity is the best perturbative quantum field theory we have 
[1] . Recently, the quantum gravitational corrections to the bending of light by the sun have been 
computed [18] and the first corrections to the black holes metrics are also known naEoj. 

It is probable that if we will ever observe quantum gravitational effects, these will be those 
already contained in the EFT of gravity, despite our knowledge or ignorance about the underlying 
fundamental theory. Since these EFT corrections are well defined, finite and computable with 
available techniques, or with foreseeable development of thereof, it becomes extremely important 
to fully determine them and to successively apply the resulting theory to any possible physical 
situation wherein an enhancement of such quantum effects could take place. Cosmology is prob¬ 
ably the best setting to look for such quantum imprints, but to be able to consistently analyze 
the EFT corrections on a Friedmann-Robertson-Walker (FRW) background, one first needs to 
compute the effective action on an arbitrary background and then look for consistent solutions 
to the effective equations of motion (EOM). In order to perform this task, we need to develop a 
covariant EFT of gravity. This is the primary purpose of this work, which focuses on the covariant 
formalism, the effective action and its curvature expansion. 

We will start by showing how the EFT of gravity can be constructed as a saddle point, or 
a loop expansion in inverse powers of the UV scale, i.e. the Planck mass, within the covariant 
formulation of quantum gravity based on the background field method formalism. We will then 
identify the classical theory (CT), the leading order (LO) quantum corrections and also present 
the form of the next-to-LO (NLO) corrections. We also outline some features of the next-to- 
next-to-LO (NNLO) corrections in the case of gravity coupled to matter. As mentioned earlier, 
it is probable that only the LO corrections will ever have a chance of being observed and so 
we will particularly focus on them, expressing all the pertaining details. Unfortunately, in four 
dimensions, it is not known how to completely compute the functional traces involved in the 
LO corrections, so we will proceed by presenting the application of the heat kernel expansion, 
both in its local and non-local forms, to their evaluation. At present, heat kernel methods 
represent the state of the art techniques to tackle these kind of contributions. We will use the 
local expansion to discuss carefully the UV divergencies and the related renormalization, which, 
as in any EFT, is performed order by order at the cost of an input phenomenological parameter 
for any divergent coupling. Then we will introduce the non-local expansion which can be used 
to evaluate the LO contributions to second order in a curvature expansion. We will discuss both 
the cases of a pure gravity EFT and its coupling to matter. In this paper we will restrict to LO 
corrections for the latter case for simplicity and leave for a future study the evaluation of those 
NLO and NNLO contributions leading to the covariant effective action that, when evaluated on 
a Minkowski background, gives the corrections to Newton’s potential obtained by the flat space 
methods discussed earlier. One particular reason behind our analysis being interesting is that 
it allows us to make a rational organization of various quantum gravity results that have been 
obtained over the years. For instance, in the EFT framework, UV divergencies can be turned into 
a physical result describing how the phenomenological parameters depend on the reference scale. 
Also the conformal anomaly is naturally included, as it induces the presence of a particular set of 
finite terms in the LO effective action, but since these terms only start to appear at third order 
in curvatures they will not be further discussed in this paper and are left to future analysis. 

The paper is organized as follows: in section]^ we develop the covariant version of the EFT 
of gravity for the case of pure gravity. In section]^ we include matter. We then discuss renormal¬ 
ization and local hnite terms in section |4] while in section [5] we describe the non-local hnite terms 
that can be inferred from the UV divergencies. In section]^ we expose the curvature expansion 
to second order and obtain the final form for our LO effective action. Finally in section we 
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summarize our results and outline their implications. In the appendix we shortly review the local 
and non-local heat kernel methods used in the computation of the LO effective action. 


2 EFT of gravity 


We don’t know the UV completion of gravity, but whatever degrees of freedom are associated 
with it and whichever symmetry characterizes it, at energies much lower than the Planck scale, 
the theory is broken to a vacuum invariant under diffeomorphisms. This vacuum is the back¬ 
ground physical geometry, described by the metric g^i,, and the remnant quantum fluctuations 
are typically suppressed by inverse powers of the Planck scale and can be treated in an effective 
manner. Thus we can construct an EFT of gravity by quantizing metric fluctuations around 
the background metric, 

dfj.L' ^ 9fiu + 


where the scale M, defined by 


M = 


^Planck 


( 2 ) 


VTOttG -v/IGtt ’ 

is related to the Planck mass Mpianck = 1/VG = 1.2 x 10^® GeV, the fundamental energy scale 
of gravitational interactions. In this paper we will often refer to M directly as the Planck mass. 


2.1 Bare action 

From an EFT point of view many terms can be added to the Einstein-Hilbert action, the first 
candidates being the curvature squared terms 12I1[22]. The most general four dimensional local 
action consistent with diffeomorphism invariance can be written as mm 


5efr[g] = j (Txy/g + C2,in -|- C2,2Rm + C2,3Rieim -|- C2,4ni? 

ao 92 


94 
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+ 0(^1 , 


(3) 
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where we have emphasized the number of derivatives present in the various curvature invariants. 
We have also defined the short hands Riem^ = R^uapR^^'^^ and Ric^ = R^^R^''. The first two 
orders and correspond to a definite number of curvatures, explicitly RP and Rp. At 0(9^), 
the two counts start to mix such that at a given order in the derivatives, we find different orders 
in the curvatures, for instance, the first three terms are R^ terms while the last one is R^ term. 
Later we will see that the local heat kernel expansion is an expansion in the derivatives, while 
the non-local heat kernel expansion is an expansion in the curvatures. 

Basically all dimensionfull couplings in the theory are written in Planck units where the pure 
numbers c* are phenomenological bare parameters. These parameters can be interpreted as the 
values of the dimensionless couplings at the Planck scale. Since their values are not known, we 
will later eliminate them in favour of the renormalized couplings which have to be measured in 
some experiment or observation at some characteristic scale much smaller than the Planck scale. 
The constants entering the Einstein-Hilbert action are obviously 

coM2 = 2A = (4) 

where A is the bare cosmological constant while G is the bare Newton’s constant. Since we are 
performing an expansion in the Planck mass we have to normalize ci = —1. The coefficient cq 
is somehow special from the EFT point of view since it comes to the left of the leading two 
derivatives interaction and is indeed enhanced by a power of rather than being suppressed 
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by powers of 1/M^. This does not pose a problem as its renormalized value is either zero or 
exceedingly small, as inferred from the present cosmological observations. Thus cqM^ ^ 1 and 
therefore, it can be considered as a mass term to the leading interactions. 

We will group all the terms of the same order in derivatives together and rewrite Q as 


5effb] =M2 


IehIo] + j^h[g\ + -j^h[g\ + ••• 


(5) 


where we have defined 

lEH[g]^M^Io[g] + h[g], ( 6 ) 

to denote the classical Einstein-Hilbert action. Note that now the action, apart the overall 
factor that we will use for the saddle point expansion, is written as a sum of terms of the 
form that makes explicit the energy/Planck mass expansion underlying the EFT 

approach. 

In d = 4, it is convenient, and physically meaningful, to rewrite the curvature square invariants 
by introducing the Euler density E = Riem^ — 4Ric^ + and the Weyl tensor, the square of 
which is = E + 2Ric^ — |i?^. The following relations allow us to switch between the Riemann 
basis {Riem^, Ric^, and the Weyl basis {C^,R^,E}, 


Riem^ =-E + 2C‘^ + 

We can now rewrite I 2 in the Weyl basis 


Ric^ = --E + -R2 , 

2 2 3 


(7) 


1 


1 


72 b] — ^C2,l + 2^2,2 + R + ) ^ “ ( o'^2,2 + C2,3 ) E + C2,4DR 

= C]i2R? + 0(^20“^ + ceE + cn/jDR. 


( 8 ) 


Using the conventions of higher derivative gravity, we find c^p. = |, 0^2 = ^ and ce ^ 

|21j . Later on we will also use the Ricci basis {Ric^,R^,E}, where one uses the Euler density 
to eliminate the Riemann tensor in favour of the Ricci tensor and Ricci scalar via Riem^ = 
4Ric2 - R2 + E. 

The six derivative term is composed of all the operators appearing in Table There are 
ten different invariants, the first two of them are of second order in the curvatures but contain a □ 
operator that makes them of 0{d^). The remaining ones have three curvatures, in four dimensions 
the last two are not independent, so there is only one invariant with three Riemann tensors, or 
equivalently with three Weyl tensors. This is the Goroff-Sagnotti invariant that characterizes the 
perturbative two loops UV divergencies of quantum gravity [221 [Ml IM]- Note that, as we will 
also discuss later, when one includes matter, the suppression factor of I 3 in ([^ has to be replaced 
by —>■ ^ 2 ^j^ 2 ; where m is the mass of the lightest particle that has been integrated out. 


2.2 Effective action 

The covariant construction of the EFT focuses on the effective action Tb], from which one then 
obtains the effective or quantum EOM 


^rb] 


= 0. 


(9) 


The solution to this equation is the vacuum or background geometry, around which small quantum 
fluctuations are quantized. It is also the condition satisfied by the on-shell metric. The virtue 
of the covariant formalism is that it allows for solutions of (j^ other than the flat Minkowski 
metric, such as static spherically symmetric or homogeneous isotropic metrics. In other words, to 
construct the EFT, we only assume that we are deep into the broken phase and the choice of the 
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Figure 1: Diagrammatic representation of the effective action of low energy quantum gravity, 
as in (13). While the dots of different colors are self-explanatory, the blue line indicates the 
propagator with the Einstein-Hilbert term. 


vacuum is not specihed. It indeed emerges as a solution to ([^ and the boundary conditions are 
ultimately provided by experiments. The effective action r[gf] is also the generating function of 
one-particle-irreducible (IPI) vertices. These are computed by taking functional derivatives with 
respect to the metric and/or matter fields and then going on-shell. As any scattering process or 
correlation function can be decomposed in IPI parts, the knowledge of the effective action allows 
their determination on any background geometry. 

As mentioned earlier, the basic underline assumption of EFT is that fluctuations are small 
which are then quantized in the standard way around an arbitrary background and this procedure 
is well under control even in the case of quantum gravity. What we lack is a theory of large 
quantum spacetime fluctuations that are expected to be large starting from the scales of 0{M) 
which will allow transitions between different vacua of the theory. The best way to quantize 
the theory keeping the background unspecified is the background field method, which has been 
applied to gravity since the early times. We will work in Euclidean signature and explain later on 
how to perform the continuation to Lorentzian signature. Within the background field method, 
the effective action is then defined by the following functional integral 


,-M2r[g] ^ 


'IPI 


e 


-Se«[9+^h] 


( 10 ) 


where 0 is implicitly used. In ( [Tol ) the background gauge fixing and background ghost are 
understood, i.e. they can be seen as part of the definitions of the gauge invariant measure 
As usual, the integral is only over IPI diagrams. More details on the background field method 
in quantum gravity can be found in [26]. In this setup the EET effective action for gravity 
is computed via the saddle point, or loop expansion in the small parameter 1/M^. We then 
substitute (j^ into (10) to find 

^-M'^[lEH[g+^h]+^l2[g+^h ]+---} 




'IPI 


e 


( 11 ) 


Expanding now the invariants in the exponential and then the exponential, all in powers of 1/M 
gives a series of Gaussian integrals, with quadratic action f yfg h ■ ■ h. These are performed 

with the aid of Wick’s theorem, leading to 


Pb] — iEH[g] + ^^h[g] + j^h[g ]— 

+ ( 12 ) 
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which can also be obtained by the standard loop expansion applied to the action We 

then expand the loops in 1/M^ and collect all terms of the same order to find 


rb] = IehIq] 


CT 


+ 


+ 


ip [h[9] + ^Trlog4^ib]} 




l^sb] + {Ieh[9\) + 2 loops with lEH[g\^ 


LO 


NLO 


(13) 


+ 0(a^) 


NNLO 


This is the final result for the covariant effective action within the EFT of gravity. It is a covariant 
expression that can be used on all backgrounds which at classical level reproduces GR, while the 
second line, the LO corrections, may contain all accessible quantum gravity phenomenology. The 
evaluation of the LO corrections is the main goal of this paper. A diagrammatic representation of 
this equation is given in Figure Before we proceed, a few comments about our main expression 
(13) are in order: 

• It clearly shows why GR is so successful: the first quantum corrections are suppressed by 
Planck mass 1/M^. 

• The fact that the propagator is given by the inverse of the Hessian of Einstein-Hilbert 
action guarantees unitarity. 


• By construction, it is invariant under gauge transformations of the background field, the 
numerical coefficients which multiply the generally covariant quantities contributing to it 
will nevertheless be different in different gauges, only on-shell Tb], or the IPI vertices 
derived from it, will be fully gauge invarianlj^ 

• In the background field approach, the meaning of going on-shell depends on the order 
considered. At LO, the background metric is the solution of the classical Einstein EOM, at 
NLO the background metric is the solution of the LO EOM, and so on order by ordei0 

• From this covariant expansion, one can immediately recover the standard rules of EFT 

[SIITIIIHIIH]. One finds that the general Lagrangian of order is to be used both at tree 
level and in loop diagrams; the general Lagrangian of order is to be used at tree level 

and as an insertion in loop diagrams; the renormalization program is carried out order by 
order. 


Thus the expansion (13) represents the basis for a well defined, consistent and predictive frame¬ 
work for low energy quantum gravity computations on an arbitrary background. It is applicable 
when the relevant energy scales are much smaller than the Planck scale E <C M, i.e. the scales 
which we can probe, or hope to probe, directly or indirectly in the near future, and where possible 
quantum gravitational phenomena may be hiding. 

Finally, if we were to restore h in equation (13), there will be a factor of h for each loop. 
For instance, the one loop diagrams in LO and NLO terms will be multiplied by h, the two loop 
diagrams in NLO will be multiplied by h?. But since both I 2 and Is are renormalized, their 
coefficients are also of order h and respectively. More precisely, the bare couplings are of 
quantum origin since they encode the information about quantum gravity in the UV. From this 
point of view, the LO corrections are of order (while the NLO corrections are of order ( 
and the EFT expansion is truly an expansion in inverse powers of the Planck mass rather than 
in powers of Newton’s constant 00 


^In particular, also the on-shell metric can be gauge dependent. For a clear example of this see m- 
^This is not the case if dealing with a perturbative or loop expansion of the partition function Z, or equivalently 
of the functional W = log Z, where on-shell always refers to the tree level or classical background. 

^For higher loops, the bare couplings are expanded similarly in powers of 
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2.3 What we know 


We can use equation (13) to reorganize many quantum gravity results by asking the question: 


what do we already know about the LO and NLO corrections? 

To start with, the local terms Ieh,I 2 , ■■■ encode the double information regarding the values 
of the phenomenological constants and their renormalization. Since these terms are local, their 
renormalized coefficients have to be measured by an experiment or observation, but the loops in 


each of the lines of (13) are UV divergent, so we first need to renormalize the theory, by absorbing 


these divergencies in the bare parameters. The freedom to choose the renormalization scale then 
leads directly to the RG running of the couplings of these invariants. The renormalizations of the 
operators in I eh and I 2 steaming from the LO loops have been studied since the works of [29] 
and lead to a first indication that quantum gravity was not perturbatively renomalizable. From 
the EFT point of view, these are instead positive results, since they tell us how the gravitational 
couplings corresponding to the lowest operators renormalize. In a later section we will study in 
detail the renomalization of these operators. 

The first finite part of the LO terms are what we can call the “leading logs”. They are 
directly related to the UV divergencies of the I 2 operators. These non-local corrections, of the 
form iilog ^R, are known since long (SOllSTlIHH|33l|3ll[35l|36l|371[38l|39] and have been discussed 
recently in the context of EFT in gnun]. But also leading non-analytical log R terms may 
be present im isl HU Haile]. We will review these terms and their derivation in section A 
curvature expansion of the LO corrections can be performed systematically by employing the 
non-local heat kernel expansion, as pioneered in [301 EH El 1121132]. In section]^ we will evaluate 
all curvature squared (7^^) terms present in the LO corrections. In principle also the TZ^ terms 
are computable using the results presented in [l8]|l9], but we will not discuss them here. 

Some of the LO terms can be obtained by integrating the conformal anomaly, but other than 
in two dimensions, this does not lead to the full knowledge of the effective action. The Reigert 
action [50] captures these contributions, but in d = 4 there is always an ambiguity in defining 
the conformal anomaly action and other options are available m- From the perspective of the 
curvature expansion the anomaly induced terms start to be present at order TZ^ and higher, even 
if matter is present [52]. In this paper we will focus on the TZ^ terms, so we will leave to a further 
work their systematic incorporation in the EFT framework. The phenomenology of these terms 
has been studied in great detail (for a review see |53] [54] and reference therein) but we stress here 
that in general these terms are not the only ones in the effective action, since those not induced 
by the conformal anomaly will also be present simultaneously. 


The Minkowski space EFT of gravity gives us indirect information about (13). The results of 


these computations are in principle obtainable from the vertices of the effective action evaluated 
on flat space. If we know the effective action to order TZ'^, by taking four functional derivatives 
with respect to the metric and setting we will be able to reproduce the four graviton 

amplitude originally obtained in m- Conversely, one can imagine to covariantize these results 
to infer the form of the LO TZ^ or TZ^ terms, as done for the TZ^ terms in [40] . 

All this was about the LO terms, what do we know about the NLO ones? Here the only 
known result is the famous computation of the two loop divergencies induced by I eh [23l (Ml 1^ 
that dictates the renormalization of the couplings in I 3 , while nothing is known about the finite 
part of these NLO corrections. 


2.4 Explicit form of LO corrections 

Before introducing matter in the next section, we take a moment to discuss explicitly the trace- 


log formula that characterizes the LO corrections in (13) for pure gravity, that from now on we 
will call T 2 = ^Trlog/g^[g']. The Hessian of the Einstein-Hilbert action is a differential operator 
and for a proper choice of the gauge fixing it is of Laplace type. This has very interesting 
consequences, as it allows the use of heat kernel techniques, both local and non-local. Although, 
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at present, we do not have the mathematical technology to fully compute the trace-log on an 
arbitrary background, exact results are available only on some specific spacetimes, like maximally 
symmetric spaces |42) . 

Skipping the details of the computations, which can be found elsewhere (see for example |26ji. 
one finds 

T 2 = ^Tr log(A 2 - 2A) - Tr log Agt , (14) 

where the differential operator for the spin two part is 

( 15 ) 

where □ = is the covariant Laplacian, 6 pa = 5 is the symmetric spin two 

tensor identity and where we have defined the tensor 


+ RX + + i?“/p) 

+ 25“^^p.) . (16) 


The effective “gravitational” mass is rm? = —2A. When we add matter (in terms of a scalar 
field (/>) there will be an extra term proportional to the effective potential V{(p) evaluated at the 
minimum v divided by M^, i.e. = —2K + V{v)/M‘^. Here and in the following, in absence of 
a cosmological constant and/or of an effective potential, we add an IR regularization mass g? as 
in the case of massless fields. For the spin one ghosts we instead have the following differential 
operator 

=-U5P; - RP;. (17) 

We remark here that these relations are valid within the de Bonder (harmonic) background gauge 
|26j . The commutators of covariant derivatives = [Vp, V^/] are 


(i«) 


The specification of the two differential operators (15) and ( |17[ ) completely determines the form 
of the finite part of the functional trace in (14) and thus the form of the LO quantum corrections 
in the EFT of pure gravity. Unfortunately, as mentioned earlier, we don’t know how to exactly 
compute these traces on an arbitrary background and one has to resort to approximations as 
those provided by the heat kernel expansion. 


3 EFT of gravity coupled to matter 

In physical cases matter is always present and therefore the EFT must also include the relative 
field fluctuations. This leads to, for example, the terms in the effective action responsible for the 
quantum corrections to Newton’s potential. 


3.1 Bare action 


For simplicity and illustrative purpose we only consider the case of a real scalar. The most general 
form of matter invariants with at most two derivatives can be written as 


Is [</>,£/] 



+ V{P) + F{P)R , 


(19) 
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where we impose i^(0) = F{0) = 0. We can use a field redefinition to eliminate Z{4>) which is 
set to unity in the following. The general form for the bare action when matter is present can be 
written as 


Ses[(l),g] = lEH[g] + -^Isifp^g] + + ^-^s,4[0,9] + -^hig] + 


( 20 ) 


where we have included the four derivative terms which will be generated by the loops. 

3.2 Effective action 


With respect to the pure gravity case, the evaluation of the effective action needs a couple of 
steps more since we need to deal with the field multiplet traces. In particular, we need to expand 
the multiplet trace in the trace-log formula in powers of the inverse Planck mas^ 

\ 


T = 


^Trlog 


^Trlog 


^Trlog 


jhh _i_ 1 , 

^EH MS-' 

^hh 

S 

1 jhx 
M^S 

A_TXh 

\ M^S 


/« 

( fhh 1 1 

^EH “T W‘ 

jhh 

0 




1 

^ M 


fhx 




Xh 


/J 




using the expansion log [^ + = \ogA+^A — ^^^[A ^B)‘^ + 0 From the relation 
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A-^B = 


\ 


( 


1 jhh 


jhh _i_ 1 jh 


\ 


0 — 

'J T-yy 
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fhx 


T' 

\ 


xh 




/ 


/ 


1 jhx 

rhh I 1 rhh-^R 


1 rXh 
tW J- q 


\ 


\ 


/ 


we see that the linear term traces to zero, i.e. tr(A ^B) = 0 , while the quadratic one gives 


ti{A-^Bf = 


rhh _ 

^EH M^-^S 


-i: 


hx 


it + -Y^I. 


I 1 Thh S tXX-^S 
“T 


xh 


fhx 


tXX^S jhh I 1 Thh S 

■' S IvB^S 


Thus, after the multiplet trace, the trace-log becomes 


T = ^Trlog 


1 


hh 


rhh I ^ T, 

^EH + Jpls 


1 . 


+ -Trlog/^^ 


1 


2 M 2 


Tr 


_ \ _ jhx ^ jXh 

Jhh I 1 Jhh S jXX^S 

-t 1717 "T TTTJ C J- ‘ 


+ 0 


M4 ) 


( 21 ) 


•^EH w 

where we used the fact that tr(A“^i?)^ = 0. Upon expanding the remaining 1/M^ in the func¬ 
tional traces we finally hnd the form of the effective action for a scalar-gravity EFT 


r = Ieh 


ct 

+ at 

r [h +Is+ \lA\oglfH + iTvlogl^t 

LO 

W M- 

r {^3 + Isa + - h^I^jhrls^jtls^ + 2 loops with Ieh\ 

V EH EH S ) 

NLO 

+ 0{ 

.A^) 

NNLO 


( 22 ) 


'‘We are expanding as {g^^, (ji) -> My). 
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Figure 2: Diagrammatic representation of the low energy EFT expansion of the effective action 
when matter, in the form of a real scalar field, is present. 


Note that Is is “Planck suppressed”, as expected by the fact that in GR the energy-momentum 
tensor of matter is multiplied by 1/M^ in Einstein’s equations. Thus classical matter appears 
at the same order as the LO quantum corrections (apart from a factor oi h). A diagrammatic 
representation of this equation is given in Eigurej^ 


The scalar-gravity NNLO terms, not shown in equation (22), will contain the second order 
expansion in 1/M^ of the first trace-log in equation (21) and the first order term of the last trace 
of the same equation together with the contribution from B)^. All these terms contribute 

to the corrections of Newton’s potential along with the LO and NLO terms [T]. In fact, after 
restoring h, the corrections to Newton’s potential are of order h/M"^ and, when compared to the 
LO corrections of order h, are suppressed by an additional factor of Thus there exist 

possible quantum gravitational effects which are less suppressed than the known corrections to 
Newton’s potential, but these cannot be observed on Minkowski space and are indeed present 
on a spacetime with non-zero curvature, as for example ERW or Schwarzschild. This peculiar 


observation is one main reason that motivates the analysis of the LO corrections in equation (22) 
and their possible physical implications. 


3.3 What we know 

As in the case of pure gravity, we can now summarize what is already known about equation 


(22). UV divergencies have been computed for many matter gravity combinations [23 ESI ESI EZ] • 
Although it was shown originally that, despite one loop gravity was finite on-shell, as soon as, 
matter was introduced this was not the case anymore. As in the pure gravity case, these results 
describe the RG flow of the relative couplings. 

As said earlier, the computation of the corrections to Newton’s potential have been performed 
on Minkowski space using Eeynman diagrams techniques [DEIEm HIKE] and further verified by 
modern methods |16j . These corrections can also be derived by the IPI diagrams obtained by 


the effective action (22) if all the relevant terms upto NNLO are considered. The polarization 


diagrams involved are covariantly described by the leading logs introduced in the pure gravity 
case, while the new non trivial part of the computation involves the matter-graviton vertex and 
four matter vertex in the presence of gravitons. The part of the leading logs of the matter- 
graviton vertex have now been computed on an arbitrary background [SB]- The full computation 
of all the relevant IPI diagrams is left for future work. 


3.4 Minimally and conformally coupled matter 


Equation (22) shows that matter induced LO corrections are approximated by the trace-log 


formula. When matter is evaluated in its vacuum configuration, it becomes as minimally coupled 
and/or conformally coupled, that we consider in the following without loss of generality. The 
bare action for a real scalar field, a Dirac spinor and an abelian gauge field is given by 


Im[4>, V’, c, c, ff] = g] +Ii['ip,g] + /i[A^, c, c, g ], 


(23) 
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where 


[V’, S'] 

h[A^,c, c, g] 


d'^xy/g I '^d^4)d^(j) + ^4>^R + -y4>^ 


(fx^'tp{'^ + m^)^) 

d^x^g I^ {d^A^^f + d^-c 


d^c 


(24) 


Here x is a parameter, only when x = 1 th® scalar action is conformal invariant (if the scalar 
has conformal weight one). The Dirac operator is defined using the covariant Dirac matrices 
7 ^ = ea 7 “ and the vierbein formalism. Note that on an arbitrary curved manifold the abelian 
ghosts do not decouple and cannot be discarded. In what follows, we will choose the gauge a = 1. 
The computation of the Hessian is straightforward and gives the following matter traces 

To = ^Trlog(Ao+ m^) 


Ti 

2 


Ti = 


--Trlog(Ai +m^) 

^Tr log Ai - Tr log(-n), 


where the Laplacians are 

Ao = -□ + -R 


Ai = -□ + -R 
2 4 


(25) 


(26) 


(Al)(; = -□50 + i^^ 

The commutators of the covariant derivatives = [V^, Vy] are 

(DW)“^ = iiV- (27) 

These last two relations are all we need to employ the heat kernel methods to compute the traces. 


= 0 




4 Local terms via heat kernel 


We are now ready to start evaluating the LO corrections appearing in (13) or in (22). The task 


amounts to evaluating the trace-log of the Einstein-Hilbert action (14), or of minimally coupled 


matter (22). First we will compute the local part of the trace-log formula. This will lead us to 


discuss the issue of renormalization and to review the well known UV divergencies of quantum 
gravity. Here we will see this problem from the point of view of the EFT of gravity, dictating 
how the gravitational couplings run. 

In both the gravitational and the matter cases we need to compute the trace-log of the form 


T = ^Trlog(A + m ^), 


(28) 


where A = — □ 1 + U is the Hessian of I eh, of the ghosts action, or of any of the matter parts 
contained in Im- In the case of massless fields, we regularize IR divergencies by adding a mass 


term m = g. The trace in (28) can be expanded in powers of the curvature by employing the 
local heat kernel expansion which we review in the appendix It is useful to rewrite the trace 
(1^ as 

/ -K{s)e-^^ , (29) 

Jo s 


T = -- 
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so that we can link it directly to the trace of the heat kernel i^(s) = Tre of the Laplacian 
operator A. We can now employ the local heat kernel expansion for a second order operator 


K{s) = 


(47r)<^/2 ^ 




(30) 


n=0 


where the i? 2 n are the integrated heat kernel coefficients, the first few are known and we have 
collected them in the appendix With the aid of this expansion, the trace can now be split 
into a divergent and a (local) finite part T = T^iv + ^finite) the first will be discussed in the next 


section, while the latter in section 4.3 


4.1 Regularization 


When the expansion (30) is inserted into (29) we encounter UV divergencies at the lower extrema 
of the integral. In four dimensions, the first three terms will be divergent, in particular they 
will have quartic, quadratic and logarithmic divergencies. To regularize (29) we will employ for 
comparison, both dimensional regularization and a cutoff regularization. 


Dimensional regularization In this case, even if we are interested in four dimensions, we 
keep d unspecified so that (29) is finite for non-integer dimensions. We finc 0 


T = - 


1 oo „ 

2(47^5o^'"'^o 


^S^n-dl2-X^-srrd ^ -^ 

0 2(47r)'^/2 


(31) 


The Gamma function r(x) has poles when x = 0, —1, —2, —3,... so we have UV divergencies for 
n — 2 < 0 for integer d. As we will see in a moment, these dimensional regularization poles in 
X = 0 , — 1 , — 2 ,... correspond to logarithmic, quadratic, quartic, ... divergencies in presence of a 
cutoff. Finite terms start from n = g + I- We now set (i = 4 — ein(31) and expand the Gamma 
functions for small e to find the following UV divergent contributions 


Tdiv = - 


2 ( 47 r) 


m* I ^ ^ ^ j Bo + - 1 + 7 j B 2 + 0 - 7 ) B 4 


(32) 


In dimensional regularization quartic and quadratic divergencies vanish for a massless theory 
m — 7 - 0 and only logarithm divergencies are left. For a massive theory we instead have all of 
them. At this point it is not clear what happens to the BgjBg,... terms in equation (31) in a 
massless theory since the m —>• 0 limit is (IR) divergent. To understand this limit one needs the 
full non-local heat kernel that we will introduce in section If we do not worry about this fact 
and dropping the finite renormalization constant we find the classical result of ’tHooft for the 
logarithmic divergencies of a four dimensional theory 


T - ^ ^ 

e( 47 r) 


-B 4 


(33) 


A theory is perturbatively renomalizable at one loop if the operators present in B 4 are also 
present in the bare action. In the EFT framework this is instead not a requirement and eventual 
divergencies will be renormalized by switching on the operators appearing in B 4 not already 
present in the bare action. 

®We recall the basic integral JJ” = r(t)/m^*. 
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Cutoff regularization We can instead apply a hard cutofl^Auv to (29) and work directly in 

d = A 

1 poo 

E / * s-Se-"*’. (34) 

2(4>r)2 ivA;„ 

Starting at n = 3 the integral is finite in the limit Auv oo and gives the same result as in the 
dimensional regularization case (31). Upon using the following integrals 


r — 

U/A?,v 


—sm 


At 


■uv _ ^^^^2 + 


A 


m 


A 2 

^'-uv 


/l/A 2 


/ 


UV 

ds 


A?,v - (log - 7 ) + O 


1/A?rv ^ 


log -7 + 0 




A 2 

^'-uv 


now gives the following UV divergent part 


Tdiv = - 


2(4 


:7r 


A^ 


uv _ + (A log % + !^ ( 3 _ ^ 


2 

A 2 

ilTT- 


+ f Auv -" 1^(1 -i)-rr? log^^ j B 2 + f log^^ - 7 ] 


(35) 


(36) 


Thus the heat kernel coefficients Bq,B 2 and il 4 are the coefficients of quartic, quadratic and 
logarithmic divergencies, respectively. Again, to renormalize the theory we need the bare action 
to contain the invariants present in Bq, B 2 , B 4 . 

A comparison between the two regularization schemes shows that, to switch from cutoff to 
dimensional regularization, one needs to replace log —>• ^ and set the other divergencies to 

zero Auv 0 and A^v 0- Note that the finite renormalization constants are also the same 
in both schemes. Clearly, the finite terms Bq,Bs,... are also independent of the regularization 
employed. 


4.2 Renormalization 


To renormalize the theory we absorb the UV divergencies into the bare couplings present in 
the actions Ieh,I 2 ,--- defining the renormalized couplings. LO renormalizes CT and LO, while 
NLO renormalizes NLO, LO and CT and so on. This will make the EFT finite at the cost 
of introducing a measured value for each divergent coupling, but it will also tell us how the 
renormalized couplings depend on the arbitrary renormalization scale. 

In the appendix [A| we evaluate the heat kernel coefficients i ?4 for all operators under consid¬ 
eration. Using ( |33[ ) we find, for the gravitational trace, the following UV divergent part in the 
Weyl basis 


^sldiv - 


1 1 

e (47r)2 



7 


4 


U9 

180 


19 

15^ 


\JR 


(37) 


In the Ricci basis the UV divergencies are instead 


e (47r) 


d'^Xy/g 


+ —R^ + —E - —DR 
10 ^ 60 45 15 


(38) 


which is the original result of |29l I59j . When using a cutoff regularization, or generally in the 
presence of masses, we also need Bq and B 2 , which we compute in appendix |A] The final result 

®Note that the UV scale Auv can be equated to M, consistently with the fact that we used the Planck scale to 
make the couplings dimensionless in the bare action, which indeed should be interpreted as the effective action at 
scale M. 
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0^12 


lC2 

1 

120 

1 

20 

1 

10 

7 

20 

IW- 

72 

0 

0 

1 

4 


1 

11 

31 

149 

lE 

360 

360 

180 

180 


1 

1 

1 

19 

lUR 

30 

30 

10 

15 


Table 1: Renormalization constants for the four derivative operators in I 2 in the Weyl basis. 


for the divergent part (modulo finite renormalization terms) is 

1 


^aldiv - - 


2 ( 47 r) 


d Xy/g 


Ayy — lOAuym^ + 5m^ log 


A 2 


uv 

2 


m 


A 2 

ilTT' 


R 


, 23,2 13 2 , 

+ |-yA„v + y™ log^ 

'20 4 180 15 ^ m? 


(39) 


where m? = —2A + V{v)/M‘^. Note that, both the graviton and the ghost contribute to the 
terms present in the massless limit while the ghost does not contribute to the massive terms. We 
also dropped the finite renormalization constants. As mentioned earlier, to switch to dimensional 
regularization just replace log —?■ 1 and set the other divergencies to zero. We remark that 

these off shell coefficients are still gauge dependent MM- 


We can now proceed to the explicit renormalization of the I eh and I 2 couplings. From (39), 


the explicit renormalization of the I 2 couplings, in the massless and dimensional regularization 
case is 




1 


c,- = c,-- 


li 


e ( 47 r )2 ’ 

where the coefficients 7 * are reported in Table This leads to the beta function^ 


gd^Ci{n) = 


li 


( 47 r) 


which can be integrated to give 


CiifJ.2) = Cj(/il) + 


li 1 
7TA2 

(47r)^ m 


(40) 


(41) 


(42) 


which relates the renormalized couplings at different renormalization scale^ This scale depen¬ 
dence of the phenomenological parameters is in principle an observable effect. 


But also A and G need to be renormalized whatever scheme is employed. One finds from (39) 
the following relations in the case of cutoff regularization 


ar 

A^ 

1 

/ 

Ur 

Ub 

47r 

f A^y - lOAi 

1 

1 

1 

( 23,2 

Ur 

Ub 


-jA?,v + 


A 2 

ml 


13 

+ —m log 
o 


ii-T- 


UV 




(43) 


^The one-loop beta function coefficient is minus the coefficient of T 

®It is clearly understood that the RG scale first used in equation and the IR regulator introduced in the 
previous sections to treat massless fields are different scales, even though we will use the same notation /r for both 
in order to keep the notation as simple as possible. 
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and in the case of dimensional regularization the following 

1 5 4 1 1 1 13 2 

Matter contributions to these relations can be worked out easily using the relative heat kernel 
coefficients computed in appendix The beta functions for A and G are then immediately 
extracted as minus the coefficient of the poles 





y;] = -^rri^ ■ 


G 


Svr 


If we now insert = — 2A we find 

fid^A = ^A^G 

07T 


^id^G = -‘^AG\ 

OTT 


(45) 


(46) 


which are the beta functions of the cosmological and Newton’s constants within the EFT ap¬ 
proach. 

As a final topic we discuss on shell UV divergencies, the only ones that could possibly be 
physical. For this we need to specify a classical background solution of Einstein’s equations. 
For simplicity we consider the case where no matter is present, so that the EOM becomes just 
= Ag^y. In dimensional regularization, the divergent part of the trace (39) then becomes 

1 If,, 58 


^aldiv - - 


(4vr) 


- / 


— a2 53^. 2 

-A^ H-Riem^ 

5 45 


where we have used G^ = Riem^ — |A^ to write it in the form as given in 
divergencies lead to the following beta function for the cosmological constant 

„ fA\ 1 58^2 


which under the assumption that A is constant [63] gives 

9Q 

fid^iAG) =- — {AGf, 

OTT 


(47) 

These on-shell 

(48) 

(49) 


so that the dimensionless coupling AG is asymptotically safe [63] • Integrating this equation gives 
a relation similar to (42). However other backgrounds have to be considered case by case. 


4.3 Local finite terms 


The finite physical part of the effective action is independent of the regularization employed 


Lfinite — 


1 


2(47r)'^/2 


E 

n=l+l 


1 


2n d 2n 2 ( 47 r )2 


1 


1 


m 


Be + O 


1 




(50) 


which is an expansion in inverse power of the mass. The local heat kernel coefficients are known, 
for any kind of matter, to order Bq [32|. In particular, for a y = 0 scalar held we hnd the following 
hnite part 


^finite 


2 (dvr) 2 m2 


d'^Xy/g 


—RUR + —RyyUR^^^ 
336 840 ^ 


+ 


1296 


-R^ - 


1080 


RR^yR^^^ - 


j:)V oct D/4 
- Ih,, -tX,, XL/v 

2835 ^ ^ “ 






+ 


45360 1620 


pa /9 p/A ^ p7 <5 
B fj, yR ^ sB a g 




( 51 ) 
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0 

1 

2 

RUR 



1 

560 



1 

840 

1 

168 



(1-x)^ 

1296 

1 

5184 

RRf^uR^^ 


1-X 

1080 

1 

1080 

RR^.uapR^’^^^ 


1-X 

1080 

7 

8640 

pM pi2 T^a 

-ft iz-Ti Q-ft ^ 


4 

2835 

25 

4536 



1 

945 

47 

7560 



1 

7560 

19 

7560 

p a/3 p ■yS n tJ.12 
^al3 -^7(5 


17 

45360 

29 

45360 

pa /3 p/r ^ p7 <5 

rt /r i/rt 7 a ^ 


1 

1620 

1 

648 


Table 2; Six derivatives invariants present in /s and their relative Bq heat kernel coefficients for 
spin 0 and ^ fields. 


The other massive matter cases can be worked out |42l [H3] and are given in Table As we noted 
before, the last two terms with three Riemann tensors are not independent in four dimensions. 
The inverse mass expansion in (50), when applicabl^ changes how six, or higher derivative terms 
are suppressed in the matter-gravity EFT —)• [2]. 

General expressions for the i? 2 n are unmanageable, but certain classes of invariants can be 
re-summed as we will see in section]^ These expressions, when expanded in the massless limit, 
give rise to leading logarithmic contributions to which we turn our attention to, in the following 
section wherein we also discuss their RG improvements. 


5 Leading logarithmic terms 


Up to now we have examined the local contributions to the effective action coming from the local 
expansion of the functional trace T. UV divergencies, being local, are the first of these terms, 
while the other terms lead to finite local corrections when the action is expandable in inverse 
powers of mass. When this is not possible, i.e. in the case of massless theories, the leading finite 
corrections are non-local logarithmic terms in d = diaolEIllMlIMlIMlEaEZj. Being second 
order in the curvatures, we can write these terms as 


^finite = J aR2R\og -^R + ac2 log 


+ .. 


(52) 


Both and a (^2 are calculable constants with the aid of the non-local heat kernel expansion, 
as we will show in the next section. But there is indeed a clever and physically transparent trick 
that allows the straightforward determination of these constants. For dimensional reasons we 


have introduced the reference scale /r in (52), which indirectly carries the information about the 


®Not applicable to gravity even if there is a cosmological constant since the ghosts are in any case massless. 
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renormalization needed to define the effective action. Notice that the loop in the trace inevitably 
leads to terms of the form log , which is then separated by introducing the scale /r to obtain 

the finite terms (52) and the divergent term — 21og^^^™, or equivalently -2- in dimensional 
regularization. Thus there is a direct link between the coefficients a, and the coefficients of the 
one loop beta functions of the couplings of the B? and invariants. A comparison with (40) 
then gives 


Oii = 


li 


2(47r) 


(53) 


and the coefficients of the logarithmic terms can be read off from Table without any further 
computation. In section we will check this result with an explicit computation. 


A comment on running couplings It is precisely the argument exposed in the previous 
paragraph that explains why and when the RG running with respect to the unphysical parameter 
jj., or equivalents, can instead be interpreted as a real, physical running with respect to changes 
of a physical scale, as can be the momenta of a photon used to “look” at a proton. Take the case 


of QED. Since the effective action contains a non-local term of the form (52), 


1 


QED 


dSvr^ 


log^T“^ + ... 


(54) 


then a variation with respect to the physical scale q and the a priori unphysical RG scale n are 
related by 

2q^dq2 = -fidfj,. (55) 

For this reason the running coupling is straightforwardly well dehned in QED or QGD where the 
couplings are dimensionless. When couplings are dimensionfull, the effective action will still be 


physically scale dependent, but one cannot rely on a simple relation like (55) to obtain it and a 
more careful analysis is indeed needed. 


RG improvement For dimensional reasons, further corrections to the leading logarithmic 
terms in (52) can be parametrized as i?/i^ 2 (—D///^) R and equivalently for the Weyl term, where 
the hi are functions of u = The logarithmic terms in (52) can be obtained by solving the 

following equation, 

uduhi{u) = aihi{u), (56) 

if we set hi{u) = 1 on the right hand side of it. The appearance of the factor hi{u) on the right 
hand side of (56) represents instead an RG improvement which, through its solution gives rise to 
the log-resummed form 


hi{u) = - 1 = ctjlogu + -(uj logu)^ + .... 


(57) 


Using this the RG improved version of (52) is then 

Rfinite — I d XyJ~g 


-□ 




1^2 


R( — ] R + 


-□ 






(58) 


where we note that 0 < Uj <C 1, even in presence of large numbers of matter fields, and so 
these corrections are genuinely non-local. Applications of this class of effective actions have been 
discussed in [66] . 

^°Not contained in the LO contribution T, but in the NLO, NNLO, ... contributions. 
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Non-local vs non—analytical Up to now we discussed non-local terms, but there is also 
the possibility that the argument of the logarithms is a scalar build out of curvatures, like R. 
Discarding the Weyl term, we will then have a contribution of the form 


/R 
d^x^/gR\og-^R+ .. 


(59) 


where 0^2 is again given by (53). This is a non-analytical contribution which can be found, for 


example, upon evaluating T on the sphere where the explicit knowledge of the spectrum of the 
Laplacian allows for a direct computation of the trace-log formula [32]. The phenomenology of 
these non-analytical terms is of the f{R) type and has been discussed extensively, see for example 
|131|331|351[16|. We here emphasise that such terms are present in the EFT of gravity with definite 
coefficients. In particular, the RG improvement can be applied also in this non-analytical case 
|44j . Before we close this section a final question to ask is: when is the action non-analytical 
and when is it non-local? The answer might depend on the properties of the background, in 
particular compact spaces usually give rise to non-analytical terms while non-compact spaces 
usually generate non-local terms I67|. In any case, both terms can in principle be present in the 
EFT effective action depending on the context, but probably not together. 


6 Non-local terms via heat kernel 


We can resum a subclass of terms in the series (50) for Tfinite by rearranging it and collecting all 


terms with the same numbers of curvatures. This naturally leads to a curvature expansion of the 
form 

^finite = 2 ~finitel 7^2 + Tfinitel7^3 + 0{TZ^) (60) 

where each term at a given order can be written in terms of the relative non-local heat kernel 


structure functions, that we shall define in the following. For example, in equation (51) the first 


term is of the form RDR, a similar term for Rg will be of the form and therefore in 

general in the coefficient B 2 n there will be a term of the form RB1^~'^R. The series so defined 
can be resummed using the non-local heat kernel expansion and gives rise to a, generally non¬ 
local, structure function of the variable □ / m?. Equivalently for the other Ti? curvature invariant 
and similarly for the other higher order curvature terms that start to appear from Rg. 
This summation approach allows us to obtain the complete form of rfinitel 7^2 and TfinitelT^s- The 
complexity of the latter is very demanding and therefore, in what follows, we will only focus on 
the curvature squared terms. 


6.1 Curvature expansion to order TZ^ 

Using the non-local heat kernel expansion reported in the appendix we find that the finite 
part of the curvature square terms in the LO effective action is given by 


Unitek.- 2{4lTy/B S 


(R^_d/2+2 


ifii-sn) - fiiO)] e 


n , (61) 


where the fi for i = {Ric, R, RU, U,Q} are the non-local heat kernel structure functions given in 


equation (90). The subtraction /i(0) reflect exactly the renormalization of the curvature square 


terms we performed in section and makes the integrals in the above equation finite in the limit 
Auv —^ oo. This allows us to define the finite TZ^ structure functions as 


7m ... 2 


X 

nv 


A^y^OO S 


(62) 


The great advantage of using this formalism is that the heat kernel structure functions do not 


depend on d and thus the relation (62) gives Tf,r,u p| 7^2 in arbitrary dimensions. 
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In d = 4 we find the following form 


^finite 1 7^2 


2 ( 47 r) 






-□ 




R 


-□ 




-□ 


+Riru[ — U + U7[/ — U + 




-□ 




(63) 


where the four dimensional finite Ti? structure functions are 


lRic{u) 

1r{u) 


lRu{u) 

lu{u) 

lQ.{u) 


40 12 u 2 Jo ^ 


- + m -0 


log[l + U^(l - 0 ] 


23 1 

960 96n 


1 

12 


-t 

2 70 


+ 2^(2 - 0(1-e')} log [l + nC(l-0] 


dC 




log[l + uOl-0] 


dC log[l + uOl-0] 


12 


d^ 




log[l + uOl-0] > 


(64) 


and u = —Hi/m?. We remark that these relations give us finite part of the effective action to order 
TZ^ for every theory whose action’s Hessian is a Laplacian of the form A = — □! + U. From these 
expressions for the structure functions one can understand the behavior of the effective action in 
the two opposite limits, namely the decoupling limit u <C 1 and the massless limit u 3> 1 . 


Local expansion 

a Taylor expansion 


In order to obtain the structure functions in the decoupling limit, we perform 
of (64) around u = 0 to find 


iRiciu) 

1r{u) 

lRu{u) 

lu{u) 

ln{u) 


u u 

15120 ~ 
u llu^ 

30240 ~ 

2 2 
U U U 


U 


30 


280 1890 


166320 

19m3 

332640 

+ 0 {u^) 


+ 0{u^ 


+ 0 (u^ 


U 4 


u u 
1680 


u 


15120 


+ O(u^) 


(65) 


Note that the linear terms in u of 'yRic and 7 r correctly reproduce the first two terms 
in equation (51). Similarly, in the case of a scalar with arbitrary x or a spinor, 
consistently obtain the results of the first two lines of Table The higher order terms 
lead to the coefficients of the corresponding operators in Hg, Hio etc. 


of order 
one can 
in u will 
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Non-local expansion In the massless limit we can make a Taylor expansion of the structure 
functions around u = oo and obtaiiF^ 


7Ric{u) 

iRiu) 

lRu{u) 

lu{u) 

lQ.{u) 


23 1 5 logM 

-log u H- 

450 60 18u Qu 


+ 


1 

1 


logtt 


2 v? 
log u 


+ 0 


_f_ 2 logu _ 

1800 120 9m 12m 8m2 4m2 


. o + 0{ ^ 


M- 


1 


1 


1 


“TT + R log^ + - - TT^ - 
18 6 M 2m^ 

1 1 log M 

1 - - log M- 

2 u u 


logM 


M^ 


+ 0( ^ 

M'^ 


1 


2 1 1 

-log M H- 

9 12 ^ 2m 


2 v? 
logM 3 
2m Av? 


+ ■5^ + 0 


log M 
2 m 2 


+ 0( ^ 


( 66 ) 


This is indeed the expansion which gives rise to non-local terms in the effective action. Note that 
the scheme dependent constant terms can be removed by a finite renormalization and thus we 


will drop them in the following. From ( 66 ) we see that in the strict massless limit, i.e. keeping 


only the logarithms, the effective action becomes 


^finite 1 7^2 — 


1 


2(47r)2 


d^xy^tr 


+ -^Rlog^R 

60 p"' 120 


-J/Jlog^U + (^Ulog^U + 

b 2 /i^ 12 




log 


(67) 


It is evident from this action that these logarithms are indeed the ones predicted in section To 
be able to check equation (53) we first need to write (67) in the Weyl basis. 


Weyl basis In d = 4 the shift to the Weyl basis is defined by 


lc{u) = ipRic{u) 

iRbisiu) = ^iRiciu) + jr{u) , (68) 


while the other structure functions remain unchanged. The shift between two basis is made using 
the generalized Euler identity [42( I67j , which also shows that the difference between the structure 
functions in different basis is of order TZ^. In the Weyl basis the heat kernel structure functions 
depend on d since both the definitions of the Weyl tensor and the Euler invariant contain it while 
in the Ricci basis they do not. Eor this reason up to now we have been employing the Ricci basis. 

The local and non-local expansions for yc and ^Rbis are 


and 


7c(m) = 

’^Rbisiu) — 


M 

1680 

17m 


+ 


+ 


M 

30240 

M^ 


M 

332640 
59m3 


5040 2592 997920 


1 


+ O(m^) 
+ O(m^) 


23 1 5 log M 

= 900 -l2o‘°®“+^- m7 + 8u2 

, , 19 1 7 logM 5 

iRbisiu) - — + — + ^ 


logM 
12 m2 


+ 


+ o{ ^ 


(69) 


(70) 


While the structure functions (641 are well defined also for negative values of u and thus of the expansion 
around u = oo imposes the restriction m? > 0 which we will implicitly assume every time we perform the massless 
limit u ^ oo. 
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In the Weyl basis, the effective action (|67|) in the massless limit therefore becomes 
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^finite 1 7^2 — 


2(47r) 


uxy/g tr 


Wlog 




( 71 ) 


We can now explicitly check the validity of (53) and therefore of the arguments presented in 
section [5l 


6.2 Effective action in = 2 

It is useful to check our relations in a case where the explicit form of the LO effective action is 
exactly known: the Polyakov action. In d = 2 we have R^y = ^g^yR, or equivalently C^yaf^ = 0, 
and thus there is only one gravitational structure function 


lR2d{u) = ]pRic{u) + 1 r{u) . 

The curvature square part of the effective action then is 
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rfinite| 7^2 — 
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(72) 


(73) 


For the case of a minimally coupled massive scalar, employing equation (62), we find 
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with the following expansions around n = 0 and u = oo 
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Similar expressions can be easily obtained for the fermions and gauge fields. Gravitons do not 
propagate in d = 2, so there is no contribution from the Einstein-Hilbert action. In the massless 
limit M —)■ oo, the last equation gives 7 R 2 d(c>o) = and we correctly recover the Polyakov action 

ESI 


If 1 

TfinitelT^a = ^ dd. 


(76) 


Note that only in the massless limit, the TZ^ or the higher order terms vanish. On the other hand, 
the massive version of the Polyakov action to order TZ^ is obtained by combining equation (73) 


together with equation (74). 


6.3 LO effective action to order 71“^ 

We can finally combine our findings for the local and non-local terms and write down the effective 
action to LO in the Weyl basis as 


ng] = 
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167rG 


d^x^{2A -R) + ^ I d^x^C^ + ^ I d^x^R^ 


+ j + j R + 0(TZ^). (77) 
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Here G,A,^ and A are the renormalized couplings or phenomenological parameters, i.e. those 
measured in experiments or observations. Thus the first line of (77) is the input to the EFT of 
gravity while the second line instead represents the output, i.e. the universal prediction common 
to all possible UV completions. Note that in (77) in order to write down the Einstein-Hilbert 
part in the conventional form we have rescaled the effective action as T —?• lOvrCT. The final 
non-local structure functions G and T are, in the cases of scalars, fermions, abelian gauge helds, 
given by 
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while in the case of gravity they are 
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(78) 


(79) 


where the first line in this last equation is the graviton contribution with = —2A + V{v)/M‘^ 
while the second line corresponds to the ghost contribution with m? = 


To make explicit the non-local terms contained in the effective action (77) we expand (78) 
and (79) in the massless limit around u = oo. For example in the case of a minimally coupled 
scalar (y = 0 ) we find 
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(80) 


while the coefficients for the other cases are reported in Table Note that the conformal invariant 
matter has vanishing leading logarithms in the Ricci scalar sector as expected. Moreover, the 
parameter y also drops out of the Weyl tensor sector. In the last grey shaded column of Table 
only the leading logarithms receive contributions from both the graviton and the ghosts, while 
the other terms, i.e. those which vanish in the massless limit, get only graviton contributions. It 
is also understood that the mass in the leading logarithms can be, modulo finite renormalizations, 
either —2A + V{v)/M‘^ or The mass in the other logarithms and in the overall mass terms 
are to be interpreted as — 2A + E(u)/M 2. In Figure]^ we have plotted the structure functions Q 
and T" of a scalar together with their small and large u approximations. From this example we 
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Table 3: Non-local coefficients for the various spins. All numbers should be multiplied by 


clearly see that the knowledge of the two expansions allows us to determine the threshold value 
of u which separates the decoupling from the massless limit, which is a priori not obvious. It also 
guides us in determining the terms which are relevant in a given regime, especially when solving 
the effective EOM. 


The action (77) is derived in the Euclidean signature while in order to understand the imprints 


in a physical context, we need to switch it to the Lorentzian signature. The final action is then 
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(81) 


where evidently the operator □ is constructed with the Lorentzian metric and more importantly, 
the Green’s function is to be interpreted as the retarded Green’s function to preserve causality 


[SnilM!. To conclude, the action ( |81[ ) is the effective action including the LO quantum gravita¬ 
tional corrections upto the second order in the curvatures, we should do physics with and it can 
be applied to any arbitrary background. 


7 Discussion and Outlook 

In this paper, we have developed the EFT of gravity in a completely covariant formalism as 
a saddle point, or a loop expansion in the inverse powers of the Planck mass scale. We have 
considered both the cases of pure gravity (with the cosmological constant) as well as gravity 
coupled to matter. By doing so, we have identified the classical theory together with the LO and 
NLO corrections and also outlined the structure of the NNLO terms in the effective action. We 
particularly focussed our attention on the local and non-local correction terms at the LO which 
were computed using the heat kernel methods. In this process, we also discussed some subtleties 
associated with the regularization and renormalization of the effective action and hnally computed 
the finite part of the effective action expanded to the second order in the curvatures. We found 
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Figure 3: Non-local structure functions (black) Q and T for a conformal scalar together with 
their small u (red) and large u (blue) approximations. The dashed lines represent the leading 
contributions while the solid lines indicate the contributions of the first four terms in each case. 
Note that the small u expansion is much less sensitive to subleading corrections than the large u 
counterpart. 


that the LO local terms in our finite effective action, consist of the and operators with 
their coefficients being phenomenological parameters fixed by observations while in the non-local 
sector, the various terms, to the second order in the curvatures, appear with their respective 
structure functions which are completely determined by the covariant EFT of gravity. 

We would like to stress that since the final effective action in our EFT of gravity has been 
computed in a complete covariant manner, it can be readily used on an arbitrary background and 
therefore, can be employed to understand different phenomenology. In particular, we have shown 
that there exist LO quantum gravitational effects which are less suppressed than the well known 
corrections to Newton’s potential but these effects are only present on a spacetime with non-zero 
curvature and cannot be observed on Minkowski spacetime. Indeed, the quantum corrections to 
Newton’s potential are of order h/M"^ and, when compared to the LO quantum corrections of 
order h, are suppressed by an additional factor of 1/M^. This fact is one main reason to study 
and understand in depth the physical implications of the LO corrections and is also the primary 
reason behind developing the covariant EFT of gravity. 

Our formalism allows us to consider the inclusion of the cosmological constant in the frame¬ 
work of the EFT of gravity as Minkowski spacetime is not a solution to the Einstein EOM in the 
presence of a cosmological constant and our ability of quantize gravity in an arbitrary background 
comes as a rescue to this. We emphasize that this is indeed an important step since observations 
indicate the presence of a non-zero cosmological constant and, even in the absence of matter, its 
presence induces all the non-local terms in the effective action which, in principle, could have a 
physically observable effect. 

The non-local terms that we have shown to appear in the effective action have coefficients 
which are completely determined by the EFT of gravity and thus, they are indeed a distinct 
prediction of our formalism. Some of these non-local terms have been recently considered in the 
literature as consistent non-local modifications of GR and have been studied in the cosmological 
context as for example candidates of dark energy in order to explain the current acceleration of 
the universe [MlEOllIIlinilTaiTaiTaiTgiTTllTSllTSlI^ in the context of black hole solutions 
|80l IHTl [82] and also in the Newtonian limit |83l 184] . Non-local modihcations of GR appear also 
in the context of super-renormalizable theories j85L ISH] . So far, these terms have been studied 
on their own together with the Einstein-Hilbert action but we have clearly shown that all such 
non-local terms appear together in the EET of gravity and therefore, should all be included in 
order to construct a viable scenario of the universe. Furthermore, we have also understood that 
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non-local terms are present in the massless limit, where the effective mass is much smaller than 
the energy scale of the process. On the contrary, in the decoupling limit, only local terms are 
present in the effective action. 

To gauge the magnitude of LO corrections in a curved background, we will specialize the 
effective action to the FRW spacetime, derive the effective Friedmann EOM and study the corre¬ 
sponding solutions, in a follow up paper lEZj. In future, we also plan to study black hole solutions 
for the effective action of the EFT of gravity and successively understand the role of conformal 
anomalies that are not considered in this paper as we restrict our analysis to curvature squared 
terms. We also intend to derive the covariant form of the effective action corresponding to the 
matter-gravity EFT, in order to obtain the known corrections to Newton’s potential. 

To conclude, the EFT formalism is a consistent approach which is able to access the low 
energy phenomenology of quantum gravity but it does not solve any of the fundamental problems 
associated with it, as it breaks down before the Planck scale. In any case, a UV completion of 
gravity capable of resolving these issues will lead to the same low energy phenomenology as the 
EFT and it may very well be, that the first quantum gravity imprints that we will ever observe 
are already described by the EFT formalism. 
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A Heat kernel methods 

When working on an arbitrary background one needs a covariant way to compute the loop di¬ 
agrams. A powerful set of techniques to perform this task is based on the heat kernel and its 
expansion |69[ H2l I88| . When dealing with one loop diagrams, as those entering the LO correc¬ 
tions, one has to evaluate only the heat kernel trace, fact that facilitates computations. The heat 
kernel Kxy{s) for a Laplace type differential operator A = —□ 1 -|- U is defined by 

{ds + ^x)Kxy{s) = 0 Kxy(0) = 5xy ■ (82) 

Formally we can write the solution as Kxy{s) = e~^^^6xy and define its trace as K{s) = Tre“^'^. 
One can thus write the trace-log as a parametric integral in s over the heat kernel trace directly 
converting any expansion of the heat kernel in an expansion for the one loop diagram. 

Unfortunately, no exact expressions for the heat kernel or its trace are known on an arbitrary 
background and so we have to resort to asymptotic expansions, which in any case are very useful. 
There are two possible expansions for the heat kernel, both of them have been used in this paper: 
the local expansion, which is an expansion in the derivatives and the non-local expansion, which 
is an expansion in the (generalized) curvatures, which we review in the following. 

A.l Local heat kernel expansion 

The local expansion of the heat kernel trace is 

. OO 

( 4 ^ ( 83 ) 

^ ^ n=0 

where the integrated heat kernel coefficients B 2 n{^) are related to the unintegrated coefficients 
52n(A) by 

R2n(A) = y d'^xtr 62 n(A) . (84) 
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Table 4: Relevant traces needed for the computation of the heat kernel coefficients Bq,B 2 and 
i?4 in d = 4. These traces also offer the expressions needed to evaluate the non-local heat kernel 
trace ^7^2 (s). 


The most useful fact about the unintegrated heat kernel coefficients is that they do not depend 
on the number of dimensions. The first three unintegrated coefficients are 


60 (A) = 1 

62(A) = if-u 
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64(A) = V - + Jnu + 

2 b b 12 

+ —Rlur.f 1 — —RIu + —R^ - —DR I , (85) 

VI8O 180 ^'" 72 30 y ’ ^ 


where = [V^, Vj^]. The next, 66(A) starts already to be quite cumbersome |32l |88]; for a 
recent computation see [65] . The coefficient 68(A) is known only in special cases. 

Using the traces reported in Table|^we can evaluate explicitly the heat kernel coefficients (85) 
in d = 4 for the various Laplacians considered in the paper. By means of the following relations 
linking the different basis, 


-Riem^ 

180 


1 

1^ 


Ric^ 




—Ric^ +-R^ +- E = -H-R^- E , 

60 120 180 120 72 360 


( 86 ) 


we can write an expression for R4(A) as 


R4(A) = j d^x^(7c2 + 7 r2 R^ + 7E R + TDR D-R) : (87) 

where the gamma coefficients are reported in Table Obviously, for photons and gravitons the 
combinations R4(Ai) — 2R4(—□) and R4(A2) — 2 B 4 (Agh) have been considered, while for the 
fermions a minus sign has been added. The heat kernel coefficients Ro(A) and R2(A) are easily 
obtained from Table HI 


A.2 Non-local heat kernel expansion 

In order to calculate the finite non-local parts of the effective action we need a more sophisti¬ 
cated version of the heat kernel expansion which resums an infinite number of local heat kernel 
coefficients. This expansion has been developed in [301 EH 113 HU Hg and retains the infinite 
number of heat kernel coefficients in the form of non-local structure functions. Keeping terms 
up to second order in the (generalized) curvatures, the non-local heat kernel expansion is 

K{s) = Kt^o{s) + K^^{s) + Kt^ 2 {s) + 0{n^) (88) 
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where K-jio{s) and K'ji{s) are just the n = 0,1 terms of the local expansion (83), while 


Kn<s) = j + lRfR{-sU)R 


+ RfRu{-sa)v + ufui-sa)\j+ . (89) 


The non-local heat kernel structure functions entering this expression are 
fmcix) = ^ + 

fnuix) = -^/(x) - ^ [/(x) - 1 ] 
fu{x) = ^/(x) 
h{x) = [/(x) - 1 ] , 


(90) 


where the basic non-local heat kernel structure function /(x) appearing in these relations is 
defined in terms of the following parametric integral 


/(x)= / 


(91) 


A simple derivation of the structure functions (90) using Feynman diagrams has been given in [67] . 


The Ricci basis non-local heat kernel structure functions (90) are independent of the number of 


spacetime dimensions, as are the local heat kernel coefficients of which they are the resummation. 
They can be used to compute all one loop bubble diagrams in any dimension for any theory 
in which the Hessian is a Laplace type operator, and probably represent the most efficient way 


to perform such computations. If we Taylor expand the non-local structure functions (90) and 


insert them in (89), they will reproduce the local heat kernel coefficients of the invariants with 


two curvatures, modulo total derivatives (SUET). 
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